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Abstract 

The phenomenology of gluon saturation at small parton momentum 
fraction, Bjorken-x, in the proton and in the nucleus is introduced. The 
experimentally-accessible kinematic domains at the nucleus-nucleus 
q ' colliders RHIC and LHC are discussed. Finally, the saturation hints 

CN ■ emerging from measurements at RHIC and the perspectives for LHC 

O ! are described. 

<D . 

1 Introduction: small-x gluons in the proton and in the nucleus 

t— I \ In the collinear factorization approach of perturbative QCD, the parton distribution functions 

(PDFs) of the proton are determined through global fits obtained using the DGLAP scale evo- 
• lution equations [1-3]. The HERA ep deep inelastic scattering (DIS) data on the proton struc- 

ture function F2(x,Q 2 ) as a function of the parton momentum fraction Bjorken-x and of the 
Qh| squared momentum transfer Q 2 , and, especially, the Q 2 slope, dF2 (x, Q 2 )/<9 In Q 2 , in the small- 

x, 3 x 1CT 5 <i<5x 1CT 3 , and smail-Q 2 region, 1.5 ^ Q 2 £ 10 GeV 2 , set rather stringent 
'— constraints on the small-x gluon distribution xg(x, Q 2 ). In this kinematic region, the gluon dis- 
tribution exhibits a strong rise towards low x and the agreement of the global fits with the HERA 
F2(x, Q 2 ) data is not as good as it is at larger values of x and Q 2 [4]. In particular, the gluon 
density xg tends to rise faster than what suggested by the data. This is due to the fact that the 
kernels of the DGLAP equations only describe splitting of one parton into two or more, so that 
the resulting evolution is linear in the PDFs. At low Q 2 , the small-x gluon density may increase 
£NJ ■ to the point where gluon fusion, gg — > g, becomes significant. Within the DGLAP framework, 

this phenomenology can be accounted for in an effective way by including nonlinear corrections 
in the evolution equations, that is, negative terms of order 0(g 2 ), 0(g 3 ), etc... that tame the 
evolution towards small x. The first nonlinear corrections, the GLRMQ terms, were derived in 
Ref. [5,6]. A more accurate description of the small-x nonlinearities is achieved in the framework 
of k t -factorization, in which the BK equation [7, 8] is used to evolve the PDFs as a function of 
x for fixed transverse momentum squared, k 2 , of the gluon. Both approaches to nonlinear gluon 
dynamics, in DGLAP and in BK, suggest that one can expected potentially-measurable effects in 
pp collisions at LHC energy, for example in heavy-flavour production [9]. 

In the case of proton-nucleus and nucleus-nucleus collisions, where nuclei with large 
mass number A are involved, the nonlinear effects are enhanced by the larger density of gluons 
per unit transverse area of the colliding nuclei. The high density of gluons at small x and small 
Q 2 induces a suppression of the observed hard scattering yields with respect to expectations 
based on a scaling with the number of binary nucleon-nucleon collisions. This reduction affects 
the kinematic region dominated by smail-rr gluons: low transverse momentum p t and forward 
rapidity y, since, at leading order, we have x ~ p t exp(— y) / ^/snn- The effect, indicated as 
nuclear shadowing, is usually accounted for in terms of a modification of the parton distribution 



functions of the nucleon in the nucleus, ff^(x,Q 2 ), with respect to those of the free nucleon, 

R A ( q 2n = Si (X,Q 2 ) 

KAX ' Q) f?(x,Q 2 ) W 

where i = q v , q sea , g for valence quarks, sea quarks, and gluons. We have shadowing, < 1, 
for x ;$ 5 x 10~ 2 . However, as we will discuss in the following, the strength of the reduction is 
constrained by existing experimental data only for x > 10~ 3 . 

The use of nuclear-modified parton distribution functions allows high-density effects at 
small x to be accounted for within the framework of perturbative QCD collinear factorization. 
However, factorization is expected to break down when the gluon phase-space becomes satu- 
rated. In these conditions, in the collision with an incoming projectile parton, the partons in the 
target nuclear wave function at small x would act coherently, not independently as assumed with 
factorization. In the limit, they may form a Colour Glass Condensate (CGC) [10]: a system, that 
can be describe in analogy to a spin glass, where gluons (colour charges) have large occupation 
number, as in a condensate. The relevant parameter in the CGC is the so-called saturation scale 
Qg, defined as the scale at which the transverse area of the nucleus is completely saturated and 
gluons start to overlap. This happens when the number of gluons, ~ Axg(x, Qg), multiplied by 
the typical gluon size, ~ l/Qg, * s ec l ua l to tne transverse area, ~ ttR\. Thus: 

2 A-xg(x,Qg) Axg(x,Qg) -j^ _ A A i/3/ / \A \ v x _ n 

Qs r T7i r A o/Q ~ A x ~ A x (Vsnn) e y , with Aw 0.3. (2) 

isR\ A z /' 3 

Q\ grows at forward rapidity, at high c.m.s. energy , and it is enhanced by a factor about 6 
(200 1 / 3 ) in the Au or Pb nucleus, with respect to the proton. Saturation affects the processes in 
the region Q 2 < Q|, where gluon recombination dominates and factorization may start to become 
invalid. 



2 Exploring the saturation region 

Figured] elaborated from Ref. [11], shows the experimental acceptances in the plane (x,Q 2 ) 
for: the nuclear DIS (lepton-nucleus) experiments NMC, SLAC-E139, FNAL-E665, EMC; 
the nuclear Drell-Yan (lepton-nucleus) experiment FNAL-E772; the RHIC (dAu) experiments 
BRAHMS and PHENIX; the experiments in preparation at LHC, ALICE, ATLAS, CMS, LHCb. 

The current knowledge of the nuclear modification of the PDFs is based on the nuclear 
DIS data, reaching down to x ^ 10 -3 . As it can be seen from the figure, the LHC will give access 
to an unexplored small-x domain of QCD. There are several model extrapolations of the amount 
of nuclear shadowing in this region, with R^ h (x, Q 2 ) ranging from 0.1 to 0.8 at x ~ 10~ 4 and 
Q 2 ~ 2 GeV 2 (see e.g. Ref. [12]). 

The estimated values of the saturation scale in heavy-ion collisions at RHIC and LHC are 
reported in the figure. For a Au nucleus probed at RHIC energy, ^/snn = 200 GeV, the estimated 
saturation scale is Qg ~ 2 GeV 2 : processes that involve gluons at x < 10 -3 -10~ 2 are affected. 
For a Pb nucleus probed at LHC energy, ^/snn = 5.5 TeV, the estimated saturation scale is Qg ~ 
5 GeV 2 : processes that involve gluons at x < 10 -4 -10~ 3 are affected. The line at Q 2 = 1 GeV 2 




Fig. 1 : The kinematic regions in x and Q 2 explored by nuclear DIS and Drell-Yan experiments, by RHIC experiments, 
and by experiments in preparation at LHC. Elaborated from a compilation in Ref. [11]. 



shows the lower limit of applicability of the perturbative QCD approach. At variance from RHIC, 
where the perturbative region and the saturation region have little overlap, at the LHC it will be 
possible to explore the saturation region with perturbative probes, like heavy quarks, and cc in 
particular. This means that discrepancies between charm production measurements close to the 
threshold and perturbative predictions could signal the onset of saturation effects. We will further 
discuss this point in Section l4~2l Another very promising approach to the investigation of small- x 
effects is by measuring hard process (jets, heavy quarks, weak-interaction vector bosons) in the 
forward rapidity region (see Section l4~TT ). 

3 Hints of saturation at RHIC 

Two experimental observations in heavy-ion collisions at RHIC support the saturation predic- 
tions of a reduced parton flux in the incoming ions due to nonlinear QCD effects. On one 
hand, the measured hadron multiplicities (see e.g. Ref. [13]), dN^/dr] « 700, are significantly 
lower than the dN^/dr] rj 1000 values predicted by minijet [14] or Regge [15] models, but are 
well reproduced by CGC approaches [16]. Assuming parton-hadron duality, hadron multiplic- 
ities at mid-rapidity rise proportionally to Q| times the transverse (overlap) area [17], a feature 
that accounts naturally for the experimentally-observed factorization of yj snn~ and centrality- 
dependences in dN^/dr] (Fig. [2], left). The second possible manifestation of CGC-like effects 
in the RHIC data is the BRAHMS observation [18] of suppressed yields of semi-hard hadrons 
(p t 2-4 GeV/c) in dAu relative to pp collisions at forward rapidities (up to rj « 3.2, Fig. [2 
right). Hadron production at such small angles is sensitive to partons in the Au nucleus with 
x mm ~ p t exp(— ??)/^/snn ~ 10~ 3 [19]. The observed nuclear modification factor, R^Au ~ 
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Fig. 2: Hints of saturation at RHIC. Left: Normalized dN^/dr] as a function of c.m.s. energy and centrality (given in 
terms of the number of nucleons participating in the collision, iVp ar t) measured by PHOBOS in Au-Au [13] compared 
with saturation predictions [17]. Right: Nuclear modification factor RdAu(pt) for negative hadrons at t) = 3.2 in dAu 
at Vsii = 200 GeV: BRAHMS data [18] compared to pQCD [19, 20] and CGC [21, 22] predictions. 



0.8, cannot be reproduced by pQCD calculations [19, 20] that include the same nuclear shad- 
owing that describes the dAu data at r\ = 0, but can be described by CGC approaches that 
parametrise the Au nucleus as a saturated gluon wavef unction [21,22]. 

4 Perspectives for LHC 

4.1 Accessing the small-x region with hard processes at forward rapidity 

The four LHC experiments — i.e. the two general-purpose and high-luminosity ATLAS and 
CMS detector systems as well as the heavy-ion-dedicated ALICE and the heavy-flavour-oriented 
LHCb experiments — have all detection capabilities in the forward direction very well adapted 
for the study of low-x QCD phenomena with hard processes in collisions with proton and ion 
beams (see e.g. Ref. [23] for more details): 

• Both CMS and ATLAS feature hadronic calorimeters in the range 3< \rj\ <5 which al- 
low them to measure jet cross-sections at very forward rapidities. Both experiments fea- 
ture also zero-degree calorimeters (ZDC, \r/\ > 8.5 for neutrals), which are a basic tool for 
neutron-tagging "ultra-peripheral" Pb-Pb photoproduction interactions. CMS has an ad- 
ditional electromagnetic/hadronic calorimeter (CASTOR, 5.3 < \rj\ < 6.7) and shares the 
interaction point with the TOTEM experiment providing two extra trackers at very forward 
rapidities (Tl, 3.1 < \t]\ < 4.7, and T2, 5.5 < |r/| < 6.6) well-suited for DY measure- 
ments. 

• The ALICE forward muon spectrometer at 2.5 < rj < 4, can reconstruct J/ip and T (as 
well as Z°) in the di-muon channel, as well as statistically measure single inclusive heavy- 
quark production via semi-leptonic (muon) decays. ALICE counts also on ZDCs in both 
sides of the interaction point for forward neutron triggering of Pb-Pb photoproduction 



processes. 

• LHCb is a single-arm spectrometer covering rapidities 1.8 < rj < 4.9, with very good 
particle identification capabilities designed to accurately reconstruct charm and beauty 
hadrons. The detector is also well-suited to measure jets, QQ and Z° — > fifi production in 
the forward hemisphere. 



4.2 Probing small-x gluons with heavy quarks 

As already mentioned, at LHC it will be possible to probe the saturation region with perturbative 
probes, such as heavy quarks. The x regime relevant for charm production in heavy-ion collisions 
at LHC (x ^ 2m c exp(— y) / ^/snn) extends down to x ~ 10~ 4 already at central rapidity y = 
and down to x ~ 10~ 6 at forward rapidity y « 4 [24]. Charm (and beauty) production cross 
sections at small p t and forward rapidity are thus expected to be significantly affected by parton 
dynamics in the small-x region. As an example, the EKS98 parametrisation [25] of the PDFs 
nuclear modification, shown in Fig. [3] (centre) for Q 2 = 5 GeV 2 , predicts a reduction of the 
charm (beauty) cross section at NLO of about 35% (20%) in Pb-Pb at 5.5 TeV and 15% (10%) 
in pPbat8.8TeV [24]. 

The comparison of heavy-quark production in pp and pPb collisions (where final-state 
effects, such as parton energy loss, are not expected to be present) is regarded as a sensitive tool 
to probe nuclear PDFs at LHC energy. The ratio of invariant-mass spectra of dileptons from 
heavy-quark decays in pPb and pp collisions would measure the nuclear modification i?g b [26]. 
Another promising observable in this respect is the nuclear modification factor of the D meson 
p t distribution, defined as: 

i? D („ r,) ~ 1 x d2Ar pVA)/ d ^ 

V(AA)(Pt,^ - ^-y d 2 N B /dptdr] ■ (3) 

In Fig. [3] (right) we show the sensitivity of -Rp Pb to different shadowing scenarios, obtained by 
varying the modification of the PDFs in the Pb nucleus (displayed, for gluons, by the curves 
labeled 'a', 'b', 'c' and 'EKS98' in the left panel of the same figure). The ALICE experiment 
will be able to measure D meson production down to almost zero transverse momentum, at 
central rapidity [27]. As shown by the projected experimental uncertainties on the D° nuclear 
modification factor in pPb, reported in the left panel of Fig. |3l this measurement is expected to 
be sensitive to the level of nuclear shadowing at LHC. 

Charmonium production at low p t and forward rapidity is another promising probe of 
small- x gluons at LHC. All four LHC experiments are expected to have good capability for J/tj; 
reconstruction in the central and in the forward rapidity region. In particular, ALICE will provide 
a measurement via di-muons in 2.5 < y < 4 down to p t [27], which probes the poorly- 
known region x < 10 -5 where current PDF parametrisations have large uncertainties. Figure [4] 
shows the J ftp rapity-differential cross section at NLO from the Color Evaporation Model [30], 
in 2.5 < y < 4 with different PDF sets, compared to the projected precision of the ALICE 
measurement [29]. 




Fig. 3: Charm production in p-Pb at LHC. Left: EKS98 [25] parametrisation for the modification of the gluon PDF 
in a Pb nucleus at Q 2 = 5 GeV 2 ~ 4mj?, along with three other different scenarios. Right: corresponding i?° A in 
p-Pb at v / snn" = 8.8 TeV and expected sensitivity of the ALICE experiment with the D° — > K~7r + measurement at 
central rapidity (|rj| < 0.9), with one year of data taking at nominal LHC luminosity [28]. 




Fig. 4: J /if) production in pp at LHC. Cross section as a function of rapidity as predicted using different PDF sets 
(details in the text) , compared to the projected precision of the measurement of the ALICE experiment [29]. 



5 Summary 

We have discussed how gluons nonlinear evolution and the phenomenology of saturation are 
expected to set on at small x in the hadrons, and how these effects are enhanced by the higher 
transverse gluon density in large nuclei. The study of this almost unexplored regime can pro- 
vide fundamental insight on the high-energy limit of QCD. We have described the experimental 
indications of the onset of saturation in heavy-ion collisions at RHIC. The LHC, as a heavy-ion 
collider, will be a unique laboratory for the investigation of the saturation regime with perturba- 



tive probes, such as forward rapidity hard processes and heavy quarks at low momentum and/or 
forward rapidity. 
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